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THERHOVISCOELASTIC  PROPERTIES  OF 
UNIDIRECTIONAL  FIBER  COMPOSITES 


1 .  INTRODUCTION 


The  present  work  is  concerned  with  evaluation  of  the  time 
dependent  thermo-mechanical  response  of  a  unidirectional  fiber 
composite  consisting  of  transversely  isotropic  elastic  fibers 
and  of  linear  viscoelastic  isotropic  matrix.  The  problem  is  of 
fundamental  importance  for  fiber  composite  structures  which  undergo 
severe  temperature  variation.  A  pertinent  example  is  a  graphite/ 
polymer  space  structure  revolving  around  the  earth  and  thus,  sub¬ 
jected  to  a  temperature  cycle  whose  minimum  is  116°K  in  the  earths 
shadow  and  maximum  589 °K  in  the  sunlight.  Under  such  conditions 
there  are  significant  variations  and  time  dependences  of  stiff¬ 
nesses,  compliances  and  thermal  expansions  and  contractions.  The 
knowledge  of  such  properties  for  the  undirectional  material  is 
essential  for  analysis  of  laminates  made  of  unidirectionally  re¬ 
inforced  laminae. 

Most  of  the  methods  of  analysis  of  viscoelastic  properties 
of  composite  materials  have  been  concerned  with  the  isothermal 
case.  The  fundamental  method  in  this  case  has  been  developed 
in  [1,  2].  It  was  shown  that  elastic  and  viscoelastic  properties 
are  related  by  the  usual  correspondence  principle  of  quasi-static 
linear  viscoelasticity.  Thus,  an  expression  for  an  elastic  pro¬ 
perty  can  be  converted  into  the  Laplace  transform  of  the  cor¬ 
responding  viscoelastic  property  which  can  then  be  obtained  by 
Laplace  transform  inversion.  Detailed  analyses  have  been  given 
in  [3]  . 

When  the  temperature  varies  with  time,  thermoviscoelasticity 
of  the  matrix  must  be  taken  into  account.  This  case  has  been 
discussed  in  [4] .  A  common  idealization  of  linear  thermovisco¬ 
elastic  behavior  is  the  so-called  thermorheologically  simple  ma¬ 
terial.  This  implies  that  viscoelastic  responses  at  any  temper- 
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ature  can  be  obtained  from  a  master  response  at  reference  tempera¬ 
ture  by  a  constant  horizontal  shift  along  the  log  t  axis.  In  such 
a  case  the  previously  mentioned  correspondence  principle  can  still 
be  retained  in  modified  form  with  respect  to  a  reduced  time  vari¬ 
able. 

However,  the  thermorheologically  simple  material  is  of  limited 
usefulness.  The  main  problem  with  this  idealization  is  that  it 
implies  that  initial  (elastic)  response  of  the  material  is  not  tem¬ 
perature  dependent  while  the  developing  time  dependent  response 
does  depend  on  temperature.  Such  an  idealization  does  not  appear 
justified  for  polymers  with  large  temperature  variations.  Conse¬ 
quently,  the  present  study  is  concerned  with  the  case  of  thermo¬ 
rheologically  complex  viscoelastic  matrix. 

In  the  following,  methods  of  analysis  will  be  developed  to 
evaluate  the  thermoviscoelastic  response  of  such  a  unidirectional 
composite  on  the  basis  of  the  composite  cylinder  assemblage  model 
(CCA)  and  the  hexagonal  array  model.  Experience  with  elastic 
analyses  shows  that  the  predictions  of  both  models  are  numerically 
extremely  close  and  agree  well  with  experiments.  The  preference 
of  one  model  over  the  other  is  a  matter  of  computational  convenience. 
The  CCA  analysis  is  very  simple  and  in  the  present  case  requires 
only  the  solution  of  simple  one  dimensional  integral  equations.  The 
hexagonal  array  model  requires  finite  element  analysis  and  is  thus 
used  only  in  those  cases  where  the  CCA  model  is  not  applicable. 


2.  THERMOVISCOELASTIC  STRESS  STRAIN  RELATIONS  OP  MATRIX 

The  stress  strain  relations  of  a  linear  viscoelastic  material 
at  constant  temperature  may  be  written  in  the  general  form 

°ij(  t)  "  cijk£(t) ek£(o)  +  fQ+  Ci  jkA  (t-t'  dt' 

(2.1) 

eij  (t)  "  sijk£(t)oki(o)  +  /  S ijki(t“t' dt' 
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where  is  the  relaxation  moduli  tensor  and  S^^ft)  is  the 

creep  compliance  tensor. 

For  isotropic  materials  it  may  be  usually  assumed  that  the 
viscoelastic  effect  under  isotropic  stress  and  strain  is  much 
smaller  than  the  viscoelastic  effect  in  shear.  Separating  stress 
and  strain  into  isotropic  parts  o  and  e  and  into  deviatoric  parts 
s^j  and  e^j  according  to  the  scheme 


°ij  “  06ij  +  sij 
eij  "  e6ij  +  eij 


a  "  1/3  °kk 
e  “  1/3  ekk 


(2.2) 


the  stress-strain  relations  become 


o(t)  ■  3Ke  ( t) 


s..(t)  -  2G(t)  ei4(o)  +  2  f t  G  (t-t ') 

13  13  Jo  3t' 
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(2.3) 


eij(t)  -  l/2g(t)si;j  (o)  +  1/2  f  g(t-t* 

•'o 


9s  i- 

) — tl  dt' 
9t 1 


where  K  is  the  elastic  bulk  modulus,  G(t)  is  the  shear  relaxation 
modulus  and  g(t)  is  the  shear  creep  compliance.  It  is  understood 
here  and  from  now  on  that  the  hereditary  integrations  commence  at 
0+ . 

In  view  of  future  developments  it  is  noted  that  a  relaxation 
modulus  or  creep  compliance  can  be  written  in  the  form 


G ( t)  -  G(o)  H ( t)  +  AG (t) 


g(t)  »  g(o)  H ( t)  +  Ag (t)  (2.4) 

g(o)  -  l/G(o) 
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where  H(t)  is  the  Heaviside  unit  step  function.  Inserting  (2.4) 
into  (2.3)  we  have 


Sij 

6i  j 


(t) 

(t) 


2G  (o)  e^ j  (t)  +  2 


r 

*  o 


l/2g (o)  si;j(t)  +  1/2 


AG(t-t') — dt* 
at ' 


SS  i  j 

Ag(t-t') — 
at' 


dt' 


(2.5) 


The  material  properties  in  (2.2-5)  are  temperature  dependent 
and  the  question  which  arises  is  how  to  represent  such  stress- 
strain  relations  when  the  temperature  varies  with  time.  This  ques¬ 
tion  has  been  the  subject  of  much  work  and  discussion  and  refer¬ 
ences  may  be  found  in  Christensen  [4]  and  Schapery  [5] . 

The  most  well  known  and  common  representation  is  the  so  called 
thermorheologically  simple  material  (TSM) .  By  this  it  1b  implied 
that  if  a  relaxation  shear  modulus,  say,  at  constant  reference  tem¬ 
perature  <f>0  is  expressed  in  terms  of  logt,  thus 

G(t,$0)  -  F (logt)  (2.6) 


then  at  some  other  temperature  4> 

G ( t ,  <J>)  «  F  [logt  +  4/  ( <#>)  ]  (2.7) 

which  implies  that  the  function  F  is  shifted  by  an  amount  1 p  along 
the  logt  axis.  With  the  definitions 

log  h($)  »  \p($) 

-t  (2.8) 

5  (t)  =  /  h [<j>  (u)  ]  du 
Jo 

the  TSM  equivalent  of  (2.3)  and  (2.5)  become 
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(2.9) 


(t) 


(t) 


-,2G(5)eij(o)  + 
=  2G(o)ei;j(C)  + 


-  l/2g(5)a  (o) 

■  l/2g (o) si j ($) 


ft  ae.  . 

2  /  G(£  -  C) — =J-  dt' 

Jo  at* 

ft  ®eii 

2  /  AGU  -  V) — dt1 

•'o  at* 


as. 

+  1/2  /  g  (5  -  5') — dt* 

at* 

as 

— il  dt1 

at* 


x 

r 


+  1/2  /  AgU  -  C) 


(2.10) 


According  to  (2.6-7)  the  relaxation  modulus  (or  creep  compli¬ 
ance)  at  any  temperature  can  be  found  by  shifting  the  corresponding 
property  at  reference  temperature  horizontally  along  the  logt  axis 
by  amount  h(<f>)  which  is  also  often  denoted  by  l/aT(<|>)  .  The  main 
problem  with  the  TSM  representation  is  that  it  dons  not  admit  tem¬ 
perature  dependence  of  the  initial  (elastic)  modulus  or  compliance. 
Indeed  it  follows  from  (2.9)  and  (2.10)  that 


si;j(o)  «  2G(o)  ei  j  (o) 


(2.11) 


Thus  the  TSM  concept  is  not  suitable  for  materials  whose  elastic 
and  viscoelastic  properties  change  with  temperature  by  similar 
orders  of  magnitude.  This  is  in  particular  the  case  for  polymers 
at  significant  temperature  variations.  More  complicated  linear 
viscoelastic  materials  are  called  thermorheologically  complex  ma¬ 
terials  (TCM) .  The  most  natural  generalization  of  the  TSM  repre¬ 
sentation  is  to  add  a  vertical  shift  in  log  property-log  t  space, 
thus  conserving  the  concept  of  a  master  viscoelastic  time  response 
which  is  shifted  to  obtain  properties.  For  discussion  see  Schapery 
[5] .  More  general  creep  type  stress-strain  relations  have  been 
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discussed  by  Harper  and  Weitsman  [6]  and  these  will  here  be  adapted 
to  the  case  of  shear.  Let  the  shear  creep  function  at  some  refer¬ 
ence  temperature  <t>Q  be  written  as 

g(t,  4>o)  “  go  H  (t)  +  Ag(t)  (2.12) 

Then  the  strain-stress  relation  for  arbitrary  temperature  varia¬ 
tion  <j>(t)  is  written  as 

(2.13) 

2fli  j  ( t)  “  9oVo{<|,)sij  (°)  +  vl(*>  f  AgU  "  ^ ' )  c)*ET  ^  V2  ( 4*  * )  j  ( t 1 )  ]  dt  * 

where  “  <j>(t)  and  »  $(t')  and  the  strain-stress  relations  in¬ 
volves  the  four  functions  VQ,  V Vj  and  h,  the  last  through  (2.8)  . 

The  creep  function  at  some  temperature  4>  T4  <l>0  is  found  from 
(2.13)  by  introducing  B^(t)  ■  si;;j(o)H(t)  and  assuming  constant 
temperature.  This  yields 

g(t,  4>)  -  V(40  [g0H(t)  f  Ag(t)]  «  V ($) g(<J>Q,  t)  (2.15) 

All  of  the  stress-strain  relations  discussed  are  representa¬ 
tions  of  varying  complexity.  The  various  shift  functions  are  intro¬ 
duced  for  reasons  of  convenient  mathematical  approximation  and  not 
on  the  basis  of  physical  arguments.  Relations  of  type  (2.13)  can 
be  written  in  analogous  fashion  for  relaxation  stress-strain  rela¬ 
tions  but  the  connections  between  relaxation  and  creep  shift  func¬ 
tions  do  not  appear  to  be  known. 

In  the  sequel  (2.13)  and  it's  inverse  will  be  written  symbol¬ 
ically 


6ij(t)  «  1/2*  ®i j  ( t) 
■ij(t)  -  2£  ei;j(t) 


(2.16) 
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where  %  and  T  are  creep  and  relaxation  operators,  respectively. 

The  stress-strain  relation  connecting  isotropic  part  of  stress 
with  isotropic  part  of  strain  is  assumed  elastic  and  therefore  has 
the  form 

t)  -  3k  l<Mt)  ]  je (t)  -  a I<j> ( t)  ] 

where  a  is  secant  thermal  expansion  coefficient  in  reference  to  a 
free  thermal  strain-temperature  diagram. 

3.  THERMOVISCOELASTIC  STRESS  STRAIN  RELATIONS  OF  UNIDIRECTIONAL 
FIBER  COMPOSITE 

3.1  General 

It  is  assumed  that  the  Unidirectional  Fiber  Composite  (UFC) 
consists  of  thermoviscoelastic  matrix,  characterized  by  the  stress- 
strain  relations  described  above,  and  of  transversely  isotropic 
elastic  fibers,  thus  carbon  or  graphite  fibers.  This  includes  as 
a  special  case  glass  fibers  which  are  isotropic.  The  fibers  are 
randomly  placed  and  therefore  the  UFC  is  transversely  isotropic 
in  the  macroscopic  sense.  It  follows,  just  as  for  elasticity, 
that  the  UFC  has  five  independent  relaxation  moduli  and  five  inde¬ 
pendent  creep  compliances.  The  former  characterize  the  response 

of  the  UFC  to  average  strains  which  are  constant  in  time  and  the 
latter  -  the  response  to  average  stresses  constant  in  time.  For 
details  see  [3] .  However,  unlike  the  isothermal  linear  visco¬ 
elastic  case  (2.1)  ,  it  is  not  possible  in  the  case  of  present 
complex  matrix  stress-strain  relations  to  write  composite  stress- 
strain  relations  of  the  composite  in  terms  of  relaxation  moduli 
and  creep  compliances.  Nor  is  it  possible  to  use  the  powerful 
Laplace  Transform  Method  to  achieve  such  a  goal.  The  specific 
form  of  such  stress-strain  relations  must  in  the  present  case  be 
uncovered  by  analysis  of  a  model  of  the  UFC. 


For  reasons  of  simplicity  and  practical  significance  we  shall 
confine  the  analysis  to  those  properties  which  enter  into  the 
stress-strain  behavior  of  a  lamina  in  a  laminate  under  conditions 
of  plane  stress.  In  the  thermoelastic  case  such  stress-strain 
relations  are 


D*  =  -  (n*a*  +  2Jl*a*) 

D*  *  -  U*a*  +  2k*a*) 
where 

E*  -  Longitudinal  Young's  modulus 
E*  -  Transverse  Young's  modulus 
v*  -  Longitudinal  Poisson's  ratio 
G*  -  Longitudinal  shear  modulus 

n*,Jl*  -  Certain  effective  moduli  of  UFC,  see  [3] 
k*  -  Transverse  bulk  modulus  of  UFC,  [3] 
a*  -  Longitudinal  thermal  expansion  coefficient  of  UFC 

a*  -  Transverse  thermal  expansion  coefficient  of  UFC 

Suppose  the  matrix  is  thermoviscoelastic  and  that  the  UFC 
specimen  is  subjected  to  constant  average  stress  and  temperature 
at  time  t=0,  thus  these  quantities  are  given  by  oi;.H(t)  and  4>II(t) 
where  H(t)  is  the  Heaviside  step  function;  then  the  strains  are 
given  by  (3.1.1a-c)  with  replaced  by  time  and  temperature  de¬ 
pendent  creep  compliances  and  aj  replaced  by  time  and  temperature 
dependent  expansion  coefficients.  Similarly,  for  application  of 
Heaviside  strains  e^H(t)  and  temperature  4>H  < t )  the  time  dependent 
average  stresses  are  now  given  in  terms  of  time  and  temperature 


dependent  relaxation  moduli  C*j  and  D*.  The  general  case  of  aver¬ 
age  strain,  stress  or  temperature  inputs  which  are  arbitrary  func¬ 
tions  of  time  is  much  more  complex.  For  constant  temperature  the 
stress-strain  relations  of  the  composite  will  be  of  form  (2.1)  but 
for  variable  temperature  this  is  no  longer  the  case  i  only  analy¬ 
sis  of  a  model  can  uncover  the  nature  of  the  stress-strain  rela¬ 
tions  . 

In  the  following  we  shall  use  two  kinds  of  geometrical  models 
to  represent  the  UFC.  The  first  one  is  the  composite  cylinder 
assemblage  which  was  first  introduced  in  [6]  and  which  has  been 
discussed  in  greater  detail  in  [3] .  This  model  consists  of  a  col¬ 
lection  of  contiguous  composite  cylinders  in  each  of  which  the 
central  cylindrical  core  represents  a  fiber  which  is  surrounded  by 
a  concentric  matrix  shell.  All  space  is  filled  out  by  composite 
cylinders  whose  diameters  diminish  from  finite  to  infinitesimal 
sizes,  figure  2.  In  all  composite  cylinders  fiber  to  matrix  shell 
radius  ratios  are  the  same. 

The  second  model  is  a  periodic  hexagonal  array  of  identical 
circular  cylinders,  figure  3.  Both  models  are  macroscopically 
transversely  isotropic.  Elastic  analyses  have  shown  that  the  ef¬ 
fective  elastic  moduli  of  these  two  models  are  extremely  close  and 
are  in  good  agreement  with  experimental  results.  The  first  model 
is  employed  in  those  cases  when  it  admits  an  analytical  solution. 

The  second  model  is  used  in  the  remaining  cases  and  must  be  analyzed 
numerically. 


3.2  Axisymmetrlc  Deformation 

We  consider  a  homogeneous  cylindrical  specimen  of  some  trans¬ 
versely  isotropic  material,  without  as  yet  specifying  its  mechanical 
nature,  which  is  subjected  to  the  displacement  field 
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U1  3  eL<t)xl 

u2  -  eT(t)x2  (3.2.1) 

u3  -  eT(t)x3 

i 

( 

(3.2.2)  | 

-1 

and  it  follows  from  the  transverse  isotropy  that  the  only  surviving 
stresses  are  j 

) 

j 

0H  •  0L(t)  a22  ■  033  ■  0T(t)  (3.2.3) 

1 

4 

which  are  related  to  the  strains  by  the  stress-strain  law  of  the  \ 

material.  Next  we  consider  any  circular  cylinder  of  radius  b  with 
generators  in  Xi  direction  within  the  specimen,  extending  through¬ 
out  its  length.  The  center  of  this  cylinder  is  located  at  point 
x°,  x°,  figure  4.  Introducing  local  coordinates  x^,x^  located  at 
the  center  and  defined  by 

x2  "  x2  X2  x3  "  x3  +  x3  (3.2.4) 

The  displacements  (3.2.1)  become 


Then  the  only  nonvanishing  strains  are 
ell  "  eL  e22  "  C33  “  eT 


U1  "  eLXl 

U2  “  eT  X2  +  ET  X2 
u3  **  eT  X3  eT  X3 


(3.2.5) 
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The  first  parts  of  the  right  sides  of  (3.2.5  b,c)  are  rigid  body 
motions  of  the  cylinder  while  the  second  parts  when  converted  to 
cylindrical  coordinates  become 

ur  »  eTr  ue  *  0  (3.2.6) 
The  stresses  (3.2.3)  converted  to  cylindrical  coordinates  are 


(3.2.7) 


and  the  rest  of  cylindrical  stress  components  vanish.  Thus  the 
cylinder  is  in  a  state  of  axisymmetric  deformation  and  stress. 

Next  we  consider  a  composite  cylinder  of  inner  radius  a  and 
outer  radius  b.  The  inner  core  is  transversely  isotropic  fiber 
material  and  the  concentric  shell  is  isotropic  (or  transversely 
isotropic)  matrix  material.  The  cylinder  surface  is  subjected 


to  the  displacements 

(3. 2. 5-6)  . 

Thus 

u^s)  *  e^, 

ur(s)  - 

o 

■ 

m 

<D 

3 

(3.2.8) 

Solving  this  axisymmetric  problem  svbject  to  the  proper  continuity 
conditions  at  interface  r=>a  we  obtain  on  r»b  a  stress  °rr(k)“0T 
and  on  the  end  sections  an  average  stress  o11*aL.  Thus,  to  an 
external  observer  the  composite  cylinder  appears  as  some  homogen¬ 
eous  transversely  isotropic  cylinder  whose  effective  stress  strain 
law  is  defined  by  the  relations  between  eL,  eT  and  oL,  o^.  This 
stress-strain  law  is  now  assigned  to  the  homogeneous  cylindrical 
specimen  which  was  our  starting  point.  It  then  follows  that  the 
surface  displacements  and  tractions  on  the  homogeneous  circular 

cylinder  of  radius  b  are  the  same  as  those  on  the  composite  cyli- 
der  of  radius  b  and  therefore  replacement  of  the  first  cylinder 
by  the  second  does  not  perturb  the  state  of  stress  and  strain  in 
the  homogeneous  specimen. 


Obviously,  the  effective  stress-strain  law  of  the  composite 
cylinder  depends  geometrically  only  on  the  ratio  a/b.  Therefore, 
any  part  of  the  volume  of  the  cylindrical  specimen  can  be  replaced 
by  non-overlapping  composite  cylinders  with  same  a/b  without  per¬ 
turbing  the  homogeneous  state  of  stress  and  strain  in  the  remain¬ 
ing  volume.  This  remaining  volume  is  made  indefinitely  small  by 
adding  more  and  more  composite  cylinders  of  diminishing  sizes. 

Thus  in  the  limit  there  is  obtained  a  composite  cylinder  assem¬ 
blage,  figure  2,  whose  symmetric  stress-strain  relation  is  that 
of  a  single  composite  cylinder.  Accordingly,  we  now  consider  a 
composite  cylinder  composed  of  transversely  isotropic  fiber  and 
thermoviscoelastic  transversely  isotropic  matrix  which  is  in  an 
ixisymmetric  state.  The  nonzero  displacements  and  stresses  are 

ur'  V  °rr'  °06  and  °zz'  They  ma*?  be  functions  of  r,  z-Xj^  and 
t  and  will  be  assigned  superscript  1  for  matrix  and  superscript 
2  for  fiber.  The  cylinder  is  subjected  to  some  uniform  tempera¬ 
ture  variation  <J>(t)  and  to  mechanical  boundary  conditions 


uz(r,L,t)  *  eL(t)L 
orr(b,z,t)  -  oT(t) 


uz(r,o,t)  -  0 


(3.2.9) 


The  interface  conditions  at  r<-a  are 
u^X)  (a,z,t)  »u^2)(a,z,t) 

°rr)  "  <^2*(a,z,t)  (3.2.10) 

uz^(a,z,t)  «uz2*(a,z,t) 


Guided  by  the  elastic  axisymmetric  generalized  plane  strain  solu> 
tion  the  displacements  will  be  assumed  to  have  the  form 


(3.2.11) 


1  -  A1(t)r  +  B1(t)/r 


ur^  "  A2(t)r 
a*1*  -  u«>  ■  eL(t)z 


where  and  A2  are  unknown  functions  of  time.  The  associated 

strains  are 


e«  ■  Ai  •  Bi/r2 


6ee>  “  Ai  +  Ve* 


e(1>  -  s<2>  -  e, 

ZZ  ZZ  h 


e(2)  m  A 
err  a2 


J2)  -  A 

Bee  *2 


(3.2.12) 


Using  the  transversely  isotropic  thermoelastic  stress-strain  re¬ 
lations  of  the  fiber  material  the  fiber  stresses  are  given  by 


0  -  -<2>  +  i.»>  ♦  t  e<2>  +  DL* 


4V  ‘  l4V  *  <k  +  °t>4V  *  ”t-GT)Eee)  +  V 


“ee’  “  ,6m  +  +  (k+8Tle99>  +  V 


(3.2.13) 


where 


UL  “  -(noL  +  2AaT) 


dT  "  “(AaL  +  2k<V 
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(3.2.14) 


\>Wu\yr  . 

(  J  1-*  .  -  u  4  <  .  «  »  '  .  ■ 


k  is  transverse  bulk  modulus,  is  transverse  shear  modulus,  n 
and  i  are  additional  moduli  and  aL,  aT  are  longitudinal  and  trans¬ 
verse  expansion  coefficients,  respectively.  All  of  these  quant¬ 
ities  may  be  temperature  dependent.  Separating  the  cylindrical 
stresses  and  strains  according  to  the  scheme  (2.2)  and  using  (2.17) 
we  have 


arr}  "  K1(2A1  +  Eo  "  3al^  +  2£(Ai/3  **  B1/r2-eQ/3) 


°ee)  "  K1(2A1+e0-3o14>)  +  2r(A1/3+B1/r2-eQ/3) 
azl]  “  K1(2A1+eo-3a1<(>)  +  2r<-2A]/3  +  2eq/3) 


(3.2.15) 


Fiber  and  matrix  stresses  must  satisfy  the  equilibrium 
equation 


fv  grr-°ee  .  0 

Tr  r 


(3.2.16) 


This  is  obviously  true  for  fiber  stresses  and  is  easily  verified 
for  the  stresses  (3.2.15). 

Introduction  of  the  relevant  results  into  (3.2.9-10)  yields 
equations  for  the  unknowns  A^,  and  Aj.  After  some  rearrangement 
these  may  be  put  into  the  form 

2K1A1(t)  +  2(l/34v2)rA1(t)-2v2rA2(t)  - 


-  oT(t)+3K1a14>(t)  -  K1eL(t)  +  2/3I>'L(t) 
2KjLA1(t)  +  8/3FA1  (t)  -  2k2A2(t)  -  2£A2  (t)  - 

-  (3Kia1+DT)<Mt)  +  (A2-K1)eL(t)+2/3reL(t) 


(3.2.17) 


which  is  a  consequence  of  (3. 2.10a)  and  (3.2.11) 

The  important  case  of  free  thermal  expansion  is  obtained  by 
subjecting  the  composite  cylinder  without  load  to  a  temperature 
step  function  rise  A$,  relative  to  reference  temperature  Then 

°L(t)  "  °T(t)  "  0  (3.2.23) 


Solving  equations  (3.2.18),  (3,2.20)  for  this  case  the  effective 
longitudinal  secant  expansion  coefficient  a*  is  defined  by  eL 
while  the  transverse  expansion  coefficient  is  defined  by  the  radial 
surface  displacement  per  unit  radius.  Thus 


aj(t,  $0,  At)  -  [A^(t)  +  B^tJ/b2]  / At 


(3.2.24) 


Another  important  case  defined  by  the  present  formulation  is 
axial  stressing  oL(t)  or  axial  straining  eL(t)  in  isothermal  con¬ 
ditions  or  together  with  temperature  changes. 


Longitudinal  Shear 


We  consider  a  unidirectionally  reinforced  specimen  which  is 
subjected  to  the  boundary  conditions 


ux(S)  -0(t)x2  u2 (S)  -  0  (t)x1  u3(S)  ■  0  (3.3.1) 


which  represents  a  pure  shear  displacement.  Then  the  only  surviving 
average  Btrain  is  E^2(t)  “  and  if  the  specimen  is  transversely 

isotropic  the  only  surviving  average  stress  is  o12(t). 

Dually,  the  specimen  is  subjected  to  the  shearing  boundary 
conditions 


T1(S)  ■  t(t)n2 


t2(s)  -  tU)*^ 


T3(S)  «  0 


(3.3.2) 


i 


h  • 


in  which  aase  the  only  surviving  average  stress  is  o12(t)  “  x(t) 
and  the  only  surviving  average  strain  is  e12(t).  It  has  been 
shown  [3]  that  for  an  elastia  fiber  composite  such  a  problem  is 
solved  by  antiplane  strain  formulation.  The  same  formulation  is 
valid  in  the  present  case.  The  displacements  in  fiber  region  R2 
and  in  matrix  region  R^  have  the  form 


u1(x)  -  4  -  0X2 

(a) 

u2(x)  -  0x2 

(b) 

(3.3.3) 

u3(x)  ■  0 

(c) 

where 

4  m 


(x2,x3,t) 

in  R 

(2) (x2#x3,t) 

in  R 

(3.3.4) 


Then  the  only  surviving  strains  and  stresses  are 


2e12"*'2 

2e13  '  *'3 

(a) 

0.,  -  rt(1) 

12  -  *  2 

°13  -  L  *  -3(l) 

in  R^ 

(b) 

(3.3.5) 

_  m  c  *(2) 

°12  °L 

°13  -  GL  ‘-T 

in  R2 

(c) 

where  r  is  defined  by  (2.17),  GL  is  the  fiber  longitudinal  shear 


modulus  and  a  comma  denotes  partial  differentiation.  The  only 
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surviving  equilibrium  equation  is 


°12'2  +  °13'3  '  0 


Introducing  (3.3.5  b#c)  into  (3.3.6)  we  have 


(3.3.6) 


V2  *  U>  - 


v2  $  (2)  -  o 


in  R, 


in  R. 


(3.3.7) 


It  follows  from  (3.3.3)  and  (3.3.1)  that  $  must  satisfy  the  bound¬ 
ary  condition 


$  (C)  -  2$x2 


(3.3.8) 


where  C  is  the  contour  of  the  specimen  section.  At  fiber /matrix 
interface  C12  the  displacements  and  tractions  must  be  continuous. 
Since  the  only  surviving  stresses  are  o^2  and  o 13  the  only  sur¬ 
viving  traction  on  C12  is 


T1  "  °12n2  +  °13n3 


(3.3.9) 


Introducing  (3.3.5  b,c)  into  (3.3.9)  and  equating  fiber  T1  to 


matrix  we  have 


3»(1)  r 

I1  JE -  -  G. 


on  C, 


(3.3.10) 


and  displacement  continuity  in  terms  of  (3.3.3)  is  satisfied  if 


If  the  boundary  conditions  are  (3.3.2),  and  assuming  that  the 
entire  boundary  C  is  in  the  matrix,  then  (3.3.8)  is  replaced  by 
a4d) 

F  — -  “  n2  (3.3.12) 

-  3n 

This  completes  the  mathematical  formulation  for  analysis  of  any 
unidirectional  specimen  under  longitudinal  shear. 

Our  purpose  is  to  analyze  the  composite  cylinder  assemblage 
model  in  terms  of  the  above  given  formulation.  A  del  .iled  analy¬ 
sis  of  elastic  longitudinal  shear  has  been  given  in  [3]  and  the 
present  situation  is  very  similar.  A  composite  cylinder  subjected 
to  (3.3.1)  will  have  on  its  boundary  tractions  of  form  (3.3.2)  and 
vice  versa.  Thus  a  composite  cylinder  behaves  to  an  external  ob¬ 
server  as  a  homogeneous  circular  cylinder  with  some  viscoelastic 
stress-strain  relation.  If  this  stress  strain  relation  is  assigned 
to  a  homogeneous  cylindrical  specimen  under  longitudinal  shear, 
circular  cylindrical  portions  of  the  specimen  can  be  replaced  by 
composite  cylinders  without  perturbing  the  states  of  stress  and 
strain  in  the  remainder.  Thus  just  as  in  the  axisymmetric  case, 
the  stress-strain  law  of  the  assemblage  becomes  that  of  one  com¬ 
posite  cylinder. 

Consequently,  we  consider  one  composite  cylinder  subjected 
to  (3.3.1).  The  boundary  and  interface  conditions  (3.3.8), 
(3.3.10-11)  for  the  two  harmonic  functions  $^(r,0,t),  ®^(r,0,t) 
then  become 

$(1) (b,8,t)  -  2 @(t)b  as  0 

4>(1)  (a,6,t)  -  $(2)  (a,0,t)  (3.3.13) 


The  functions  which  solve  this  problem  are 


4>  (r,8,t)  ■  |A1(t)  +  B1(t)/r|  cos0 


4>^(r,8,t)  »  A2(t)  cos© 

(3.3.14) 

where  A^ ,  A2  and  B2  are  defined  through 

(3.3.13) 

by  the  equations 

h1(t)  m  ^2(3  (t)  -  VjAjtt)]^ 

(a) 

Bx(t)  -  a2  |a2  (t)  -  2(3(t)j/Vl 

(b) 

(3.3.15) 

(l+v2)  r  A2(t)  +  VjGjAjft)  -  4  T3(t) 

(a) 

The  last  equation  is  an  integral  equation  which  can  be  sim¬ 
plified  by  operating  on  both  sides  with  ^  which  is  the  inverse  of 
£1  see  (2.16).  Remembering  that  B(t)»e12(t)  we  have 

(l+v2)  A2(t)  +v1GLXA2*t*  "  4®i2^  (3.3.16) 

The  boundary  traction  T1(b,8,t)  is  the  left  side  of  (3. 3.13o) 
evaluated  at  r-b.  It  follows  from  (3,3.14)  that  it  is  given  by 

T1(b,e,t)  -  r  (a1  -  s1/b2)  cose  (3.3.17) 

If  the  cylinder  were  homogeneous  the  boundary  tractions  would 
have  the  form  (3.3.2).  Since  cos0«n2  on  r*b  it  is  seen  that  (3.3.17) 
is  of  that  form  (T2  vanishes  on  r»b  since  n^"0  but  does  exist  on 
the  end  sections  where  n2»0,  n^"l)  .  Therefore  the  average  stress 
response  o12(t)  is  identified  by  the  multiplier  of  cos8  in  (3. 3. IB). 
Using  (3.3.16)  we  have  after  rearrangement 


where  A2  is  defined  by  (3.3.16).  If  desired,  can  be  eliminated 
by  use  of  (3.3.16)  and  then  (3.3.18)  assumes  the  operational  form 

( 1+ V2+v1Gl^ ) a 1 2 ( t ) =2  [ vx r+ ( 1+ v2 >GL]e12(t>  (3.3.19) 

Equation  (3.3.19)  defines  the  stress-strain  relation  of  the  fiber 
composite  as  modelled  by  the  CCA.  The  relaxation  modulus  is 
given  from  (3.3.19)  as  o^2(t)  when  e12(t)  ■  H(t)  and  dually  the 
creep  compliance  is  given  as  e12(t)  when  o12 (i)”H(t)  .  When  the 
matrix  is  elastic  and  ^-1/G^  and  the  known  results  for  the 

elastic  longitudinal  effective  shear  modulus  [8,3]  are  obtained. 

3.4  Transverse  Shear  and  Transverse  Uniaxial  Stress 


Let  the  composite  be  subjected  to  the  average  stress  state 


000 
o  c22 (t)  o 
000 


(3.4.1) 


The  loading  which  produces  (3.4.1)  is  shown  in  figure  5.  The 
stress  (3.4.1)  may  be  split  in  the  following  fashion 


1 

0 

0 

0 

■ _ 

1 - 

0 

0 

0 

1 _ 

m 

0  ^0^2 (t)  0 

+ 

0  *i022(t)  0 

°  0 

0  0  “^°22(t). 

The  first  part  is  an  isotropic  transverse  stress  which  has  been 
considered  in  section  3.2  while  the  second  is  a  pure  transverse 


shear  of  magnitude  t (t) =*so22 oriented  at  45°  with  respect  to 
the  material  axes  x1#  x2,  x3<  Since  the  material  is  linear  the 
response  to  transverse  uniaxial  strain  is  obtained  by  superposi¬ 
tion  of  the  isotropic  stress  response  and  the  transverse  shear 
response . 

Unfortunately,  however,  the  CCA  model  cannot  be  directly 
analyzed  for  transverse  shear.  This  is  a  well  known  difficulty 
which  arises  also  in  the  elastic  context.  Therefore,  a  different 
model  is  needed  for  this  case.  A  suitable  model  is  hexagonal  ar¬ 
ray  of  identical  circular  fibers.  Such  a  model  has  the  required 
transverse  material  isotropy.  It  is  well  known  that  in  all  cases 
where  CCA  analyses  can  be  performed  hexagonal  array  and  CCA  re¬ 
sults  are  numerically  extremely  close  (see  i.e.  [3]),  therefore 
no  practical  inconsistency  is  expected  by  reverting  more  to  the 
hexagonal  array  model. 

Analysis  of  the  periodic  hexagonal  array  is  performed  using 
finite  element  analysis  techniques.  The  finite  element  developed 
for  this  analysis  is  a  four  noded  linear  rectangle  of  the  seren¬ 
dipity  type.  The  inclusion  of  a  linear  displacement  function  in 
the  direction  (independent  of  x2,x3)  yields  the  capability 
for  generalized  plane  strain  analysis.  Viscoelastic  material 
properties  are  incorporated  into  a  pre-stress  as  described  in 
the  following  section. 

The  periodic  hexagonal  array  (figure  3)  is  such  that  only  a 
single  repeating  element  need  be  modelled.  The  various  symmetrices 
of  the  model  and  loading  allow  for  the  prescription  of  the  neces¬ 
sary  boundary  conditions  on  the  repeating  element.  The  appropriate 
boundary  conditions  for  the  repeating  element  can  be  determined 
by  noting  that  each  of  the  hexagonal  cylinders  are  identical  and 
that  under  any  loading  state,  the  array  of  hexagons  must  be  con¬ 
tiguous  . 

In  the  case  of  transverse  (x2)  extensions!  or  compressional 
loading,  the  loading  and  model  symmetrices  require  that  sides  1, 

3,  and  4  of  the  repeating  element  (figure  3)  remain  straight, 
and  parallel  to  their  original  positions.  Contiguity  between 


hexagonal  cylinders  prescribes  that  along  side  2,  the  line  con¬ 
necting  the  fiber  centers  define  a  point  of  cylic  symmetry.  There 
fore,  points  on  side  2  above  the  intersection  of  the  fiber  center 
line  deflect  exactly  opposite  of  equidistant  points  on  the  other 
sides.  Stress  boundary  conditions  are  incorporated  in  the  usual 
fashion  with  the  stresses  being  integrated  to  determine  applied 
nodal  forces . 

4.  NUMERICAL  PROCEDURES  AND  IMPLEMENTATION 

The  analytical  models  developed  in  the  previous  section  in¬ 
volve  the  integral  operator  representation  of  the  matrix  shear  re¬ 
sponse.  The  presence  of  the  integral  operator  necessitates  that 
the  models  be  evaluated  incrementally  in  time.  The  procedure 
utilized  for  both  the  CCA  and  hexagonal  array  models  involves 
specifying  the  time  history  of  applied  load  and  temperature  and 
determining  the  response  of  the  composite  at  timeB  0,  At,  2At,  3At, 
etc.  where  At  is  an  increment  of  time.  At  each  time  step,  the 
hereditary  integrals  are  evaluated  numerically  using  simple  trap¬ 
ezoidal  integration. 


4.1  Composite  Cylinder  Assemblage 

Literature  review  has  yielded  matrix  thermoviscoelastic  ma¬ 
terial  properties  in  the  form  of  compliance  data  only.  The  equa¬ 
tions  3.2.17,  18  and  3.2.20  have  therefore  been  converted  from 
relaxation  to  compliance  form  utilizing  the  relation  %  £x»x.  This 
operation  yields  the  equations 

2^;  [kiAi  (■*=•)]  +  2(1/3  +  v2)A1(t)-2v2A2  (t)» 

^C1  (t)  +3](  [k^a^tji  (t)  ]-]£  (t)  ]  +2/3el  (t)  (4.1.1) 
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21[k1A1(t)]  +  S/SA^t)  -  2x[k2A2(t)]  -  2A2(t) 


I 

I 


O 


) 

A 


’•  •< 


3x  [  (k1a1+DT)  <t>  (t)  ]  +x  [  U2  (t)  -kx)  eL  (t)  ]  +2/3eL (t)  (4.1.2) 

2vl^  tklAl  3  ~^/3v^A^  (t)  +2v2^  [A2A2  (t)  ]  +vx^  [k^  (t)  )+v2^  (n2eL  (t)  ] 

+4/3v1eL(t)=^  aL(t)+3v1^[k1a14i(t)]~v2^[DL<|>(t)]  (4.1.3) 

where  both  the  matrix  bulk  modulus  k^  and  the  fiber  elastic  con¬ 
stants  may  be  temperature  dependent  and  therefore  dependent  upon 
the  time  history  of  the  applied  temperature. 

During  the  first  time  step  evaluation,  the  hereditary  integral 
operator  reduces  to  simply 

](X(o)-g(o)x(o)  (4.1.3) 

and  the  problem  reduces  to  an  elastic  analysis  with  the  initial, 
t»o,  loadings  and  the  initial  matrix  shear  compliance,  g(o).  Sub¬ 
sequent  time  steps  are  evaluated  utilizing  a  trapezoidal  numerical 
integration  scheme  where  the  hereditary  integral  operator  is  eval¬ 
uated  as 


][x«C1+C2x(ti)  (4.1.4) 

The  constants  and  C2  are  functions  of  the  prior  time  history 
of  x(t)  and  cj>(t).  This  yields  three  linear  simultaneous  equa¬ 
tions  for  A1(t^),  A2(t^),  and  Thus  the  problem  is  easily 

solved. 

4.2  Finite  Element  Analysis 

Implementation  of  the  thermo-viscoelastic  matrix  constitu¬ 
tive  relations  within  the  finite  element  analysis  is  accomplished 
in  the  same  manner  as  in  the  CCA  model.  The  viscoelastic  strain 
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(or  stress)  is  simply  treated  as  a  pre-stress  in  the  analysis. 

The  form  of  the  thermoviscoelastic  constitutive  relations 
for  the  matrix  were  defined  in  equations  2.3  with  the  shear  hered¬ 
itary  integral  defined  in  equation  2.13.  Utilizing  the  compliance 
form,  an  elemental  constitutive  matrix  is  formulated 

|e}“  [D]“1|o| 

where  [D]”1  is  viscoelastic  and  the  compliance  matrix  ID]'"1  is  com¬ 
prised  of  terms  involving  the  inverse  of  the  matrix  bulk  modulus 
operating  on  the  isotropic  part  of  the  stress  tensor  and  the 
shear  compliance  operator  ^  operating  on  the  deviatoric  part  of 
the  stress  tensor  s^ j .  Evaluation  of  the  compliance  matrix  at 
any  time  step  is  therefore  evaluated  using  the  same  simple  trape¬ 
zoidal  integration  scheme  as  in  the  CCA  model  yielding  . 

{e  (t±)}  -  IDJ1  (0(^)1  +{tp|  (4.2.1) 

t 

This  form  is  equivalent  to  equation  4.1.4  and  the  constant  matrix 
[D.^]”1  and  vector  jtp|  are  functions  of  the  prior  time  history  of 
4>(t)  And  jo(t)|  .  The  constitutive  relations  4.2.1  are  then  inverted 
to  yield 


(<7(^)1=  tDi]  je  (t±)f  +  jopJ 

The  development  from  this  point  on  is  the  standard  method  for  ob¬ 
taining  the  elemental  stiffness  and  load  vector  for  finite  element 
analysis  and  will  not  be  detailed. 

5.  PREDICTED  COMPOSITE  BEHAVIOR 

The  two  analytical  models  developed  in  the  previous  sections 
have  been  utilized  to  predict  the  response  of  a  unidirectional 
fiber  composite  with  a  thermo-viscoelastic  matrix.  The  results 
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of  the  analysis  are  presented  in  this  section. 


5.1  Material  Properties  and  Finite  Element  Mesh 

The  fiber  properties  utilized  are  typical  carbon  fiber  data 
and  are  given  in  table  1.  It  has  been  assumed  that  these  properties 
are  independent  of  temperature. 

Matrix  properties  utilized  in  the  analysis  are  given  in  table 
1  and  figures  6-9.  These  properties  are  derived  from  data  obtained 
from  [6]  utilizing  some  basic  assumptions  and  modifications.  The 
data  obtained  from  [6 } included  a  complete  characterization  of  the 
extensional  properties  of  a  3502  matrix  material  system  over  the 
temperature  range  303°K-423°K.  The  thermal  expansion  coefficient 
was  found  to  be  constant  ever  this  temperature  range.  Since  two 
material  constants  are  needed  to  define  an  isotropic  material  and 
only  extensional  data  were  given *  Poisson's  ratio  at  the  reference 
temperature  was  assumed  to  be  0.3*  thus  defining  the  matrix  initial* 
reference  temperature  properties.  The  time  dependent  portion  of 
the  matrix  shear  compliance*  Ag*  is  shown  in  figure  6  and  was  com¬ 
puted  directly  from  the  data  in  [6]  noting  that  all  time  depen¬ 
dent  extensional  response  is  due  to  shear  alone. 

Determining  the  temperature  dependence  of  the  matrix  para¬ 
meters  gQ*  h*  V ^  and  v2  required  modifications  to  the  data  from 
[6] .  The  data  were  obtained  in  a  form  which  could  not  be  used  be¬ 
low  the  reference  temperature,  303°K.  Additionally*  the  two  ver¬ 
tical  shift  functions*  v.^  and  V2  were  not  well  balanced  when  ex¬ 
trapolated  beyond  the  temperature  range  over  which  they  were  ob¬ 
tained.  These  difficulties  and  the  lack  of  a  second  matrix  material 
constant  as  described  above  were  remedied  by  modifying  the  form  of 
the  mathematical  representation  of  all  temperature  dependent  ma¬ 
terial  parameters  which  were  obtained  from  [6] .  After  modifi¬ 
cation*  a  single  vertical  shift  function  was  obtained  as  the  pro¬ 
duct  of  the  functions  V,  and  V«  and  this  function  was  used  to  re- 

J.  4* 

place  V2  in  equation  2.13  with  set  to  unity.  Finally,  it  was 
assumed  that  the  matrix  bulk  modulus  was  independent  to  tempera- 


ture.  These  modifications  yielded  the  data  shown  in  figures  7-9. 
While  these  data  are  no  longer  purely  experimentally  obtained, 
they  do  provide  a  consistent  set  of  properties  for  analysis. 

In  both  the  CCA  and  finite  element  analytical  models,  the 
fiber  volume  fraction  was  set  as  0.6.  The  finite  element  model  of 
the  repeating  element  of  the  periodic  hexagonal  array  is  shown  in 
figure  10.  The  loadings  for  the  finite  element  analytical  model 
consisted  of  transverse  extension  and  temperature.  Boundary  con¬ 
ditions  for  these  analyses  were  discussed  in  section  3.4. 

5.3  Axial  Shear  and  Temperature  Loading 


The  impetus  for  developing  the  analytical  models  described 
earlier  was  to  determine  the  response  of  a  unidirectional  composite 
subjected  to  combined  loading  and  cylic  temperature.  The  temper¬ 
ature  cycle  of  interest  has  extremes  of  116°K  and  589°K.  Before 
performing  the  analysis  corresponding  to  such  a  thermal  cycle,  iso¬ 
thermal  analysis  at  the  two  extremes  and  at  the  reference  temper¬ 
ature  were  performed. 

The  CCA  model  was  utilized  with  a  constant  axial  shear  stress 

of  1  MPa  at  temperatures  of  116°K,  303°K,  and  589°K.  The  resulting 

axial  shear  strains  are  shown  in  figures  11  and  12.  These  data 

demonstrate  the  acceleration  of  creep  strain  and  softening  of  the 

initial  and  time  dependent  portions  of  the  composite  response  as 

temperature  is  increased.  The  axial  shear  strain  at  116 °K  appears 

to  be  nearly  constant.  This  phenomenon  is  a  function  of  the  hor- 
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izontal  shift,  h,  approaching  10  and  the  vertical  shift,  V2* 
approaching  0.5  at  this  low  temperature.  The  very  small  value  of 
"h"  causes  the  creep  rate  to  decrease  so  that  in  effect,  the  creep 
strain  after  one  hour  is  comparable  to  a  reference  temperature 
creep  strain  after  lO”4  hours.  The  vertical  shift  factor  at  116 °K 
has  the  effect  of  nearly  halving  the  creep  strain  response  that 
would  occur  at  lo”4  hours  at  room  temperature.  The  combination 
of  these  two  shift  functions  cause  the  resulting  creep  strain  to 
appear  constant  in  the  figure. 


The  thermal  strains  produced  by  step  temperature  changes 
from  the  reference  temperature  to  116®K,  304°K  and  589°K  are  shown 
in  figure  13.  The  strains  have  been  normalized  with  respect  to 
the  temperature  change  yielding  secant  thermal  expansion  coeffic¬ 
ients  at  these  three  temperatures.  These  data  demonstrate  the 
same  horizontal  and  vertical  shift  effects  as  were  seen  in  the 
axial  shear  compliance.  The  strains  at  304°X  and  116°K  are  re¬ 
drawn  with  an  expanded  scale  in  figure  14.  The  data  again  show 
the  very  low  temperature  response  to  be  nearly  constant  in  time. 

Having  determined  the  composite  response  to  step  mechanical 
and  temperature  loadings ,  an  analysis  was  performed  utilizing  a 
constant  axial  shear  stress  of  1  MPa  in  conjunction  with  the  tem¬ 
perature  cycle  shown  in  figure  15.  The  24  hour  period  can  be  con¬ 
sidered  representative  of  the  temperature  cycle  experienced  by 
satellite  in  geosynchronous  orbit. 

The  axial  shear  strain  produced  by  this  loading/oyclic  tem¬ 
perature  combination  over  48  hours  (2  periods)  is  shown  in  figure 
16.  The  shear  strain  is  seen  to  increase  sharply  during  the  init¬ 
ial  8  hours  while  the  temperature  is  increasing.  As  the  temper¬ 
ature  drops,  the  strain  decreases  reflecting  the  increased  matrix 
initial  shear  modulus  at  lower  temperatures.  The  difference  in 
shear  strain  at  24  and  0  hours  is  the  accumulated  creep  strain 
after  one  temperature  cycle. 

During  the  second  temperature  cycle,  the  response  is  similar 
but  the  increase  in  creep  strain  is  much  smaller  indicating  that 
at  some  future  time,  the  composite  response  may  settle  into  purely 
periodia  steady-state  behavior. 

The  free  thermal  expansion  strains  produced  by  the  cylic  tem¬ 
perature  are  shown  in  figures  17  and  18.  The  evidence  of  visco¬ 
elastic  behavior  in  these  results  is  the  non- zero  strain  at  24 
and  48  hours  where  the  temperature  is  at  the  reference  temperature 
and  no  temperature  change  is  applied.  The  increment  of  visco¬ 
elastic  strain  induced  between  24  and  48  hours  is  seen  to  be  much 
smaller  than  between  0  and  24  hours.  As  with  the  axial  shear 
strain  response,  this  would  indicate  that  a  steady  state  response 
may  occur  after  further  cycling. 

-29- 


The  total  effects  of  the  thermoviscoelastic  matrix  material 
on  pure  thermal  expansion  is  considerably  smaller  than  for  the 
mechanical,  axial  shear  load.  This  is  due  to  the  more  nearly  hydro¬ 
static  matrix  stress  state  induced  by  a  temperature  change  alone. 

5.4  Transverse  Extension  and  Temperature  Loading 

The  finite  element  analytical  model  was  also  used  to  predict 
the  composite  response  under  combined  mechanical  and  cyclia  tem¬ 
perature  loading.  The  mechanical  loading  consisted  of  a  1  MPa 
transverse  extensional  stress  while  the  temperature  cycle  was  as 
shown  in  figure  15. 

The  composite  strain  response  to  this  loading  history  is 
shown  in  figure  19.  The  magnitude  of  the  applied  stress  is  small 
and  the  overall  composite  response  is  dominated  by  free  thermal 
expansion.  When  the  free  thermal  expansion  strains  are  removed, 
the  effects  of  the  thermoviscoelastio  matrix  on  the  composite  re¬ 
sponse  are  readily  apparent  as  seen  in  figures  20,  21,  and  22. 

These  figures  persent  the  mechanical  strain  in  loading  direction, 
(x2)  in  the  direction  perpendicular  to  the  load  (x3)  and  in  the 
axial  direction  (x1) ,  respectively. 

The  composite  response  without  thermal  expansion  is  seen 
to  reflect  the  same  phenomena  that  were  present  under  axial  shear 
and  thermal  cycling.  The  response  during  the  Becond  temperature 
cycle  is  similar  to  that  exhibited  during  the  first  cycle  but 
with  a  reduced  total  change  in  strain  over  the  period  of  the  cycle. 

6.  DISCUSSION  AND  RECOMMENDATIONS 

The  development  of  thermoviscoelastio  analytical  capabilities 
for  unidirectional  fiber  composites  has  provided  a  means  for 
realistic  modelling  of  material  behavior  under  cyclic  loading  and 
temperature  environments.  These  capabilities  have  been  demon¬ 
strated  through  the  prediction  of  composite  response  under  com¬ 
bined  Heavieside  mechanical  load  and  cyclic  thermal  load.  The 
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results  have  demonstrated  the  significance  of  temperature  depen¬ 
dence  of  the  various  material  parameters  in  the  matrix  represen¬ 
tation. 

The  predicted  composite  compliances  under  cyclic  temperature 
have  indicated  that  significant  differences  in  response  occur  during 
the  first  and  second  temperature  cycles  while  indicating  that  under 
subsequent  cycling,  the  material  response  may  become  steady-state. 
This  area  will  require  further  study  to  adequately  define  the  ma¬ 
terial  response  to  cyclic  temperature. 

The  difficulties  encountered  in  acquiring  matrix  thermovisoo- 
elastic  material  properties  and  the  sensitivity  of  the  predicted 
composite  behavior  have  demonstrated  the  need  for  further  study 
of  the  proper  form  of  the  hereditary  integral  representation.  This 
area  also  includes  the  representation  of  the  various  temperature 
dependent  parameters  within  the  hereditary  integral. 

The  analytical  models  developed  have  proven  the  significance 
of  and  need  for  the  thermorheologically  complex  matrix  represen¬ 
tation.  The  much  more  realistic  matrix  representation  utilised 
here  has  provided  a  pioneering  effort  in  the  analysis  of  fiber 
composites . 
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Table  1 


Material  Properties 


Fiber  Properties 
EA  -  231  GPA 

s 

E(J,  -  22.4  GPa 

G.  «  22.1  GPa 
A 

G^  “  8.30  GPa 
vA  -  0.30 
aA  ■  -1.33  ura/m°c 
aT  -  7.04  um/m°c 
Vp  -  0.60 

Matrix  properties  at  reference  temperature  (303°K) 

k  -  3.58  GPa 

g(o)  -  0.606  (GPa)-1 

v  ■  0.30 

a  ■  30.0  ym/m°K 
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Figure  1.  Unidirectional  Fiber  Composite 
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Figure  6.  Time  Dependence  of  Matrix  Shear  Compliance 
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Figure  11. 


Axial  Shear  Strain  for  1  MPa  Shear  Stress 
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Transverse  Thermal  Strains  Divided  by  Tem¬ 
perature  Change 
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Figure  16 .  Axial  Shear  Strain  Due  to  1  MPa  Stress  and 
Cyalic  Temperature 
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Figure  17 .  Free  Thermal  Expansion  Strain  for  Cyclic 
Temperature 
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Strain  History  for  1  MPa 
Cyclic  Temperature 
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